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action with the larger environment of cosmic medium. In this lecture we first review 
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0\ I Quantum dynamics of a finite degrees of freedom is often much affected by inter- 

p ■ some recent developments of the theory of quantum dissipation in the linear open 

(-H I system. In the second part we discuss two classes of applications: decay of unstable 

particle in medium, and environmental effect on the parametric particle production. 

• ^H , 

^ ■ The first subject of particle decay may have important consequences on the scenario 

H ■ 

Cu I of baryogenesis due to a GUT heavy particle. On the other hand, the parametric 

resonant particle production is related to the reheating problem after inflation, and 

its environmental effect is related also, perhaps more importantly, to the decay of 

moduli fields in supergravity models. 
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I Introduction 

Ideal elementary process rarely occurs, presumably except in carefully prepared 
laboratory experiments. It might be useful to recall that even processes usually con- 
sidered elementary can be regarded taking place in complex environments. This has 
to do with how one separates a system in question from a surrounding environment. 
For instance, the beta decay, the fundamental weak process oi n -^ p + e + Ue, when 
it occurs in nuclei, is compounded by effects of nuclear strong interaction with the 
rest of nucleons. Presumably the most spectacular of this kind is how nucleon decay 
proceeds. The extremely weak process of baryon number violating process at the 
quark-lepton level must inevitably occur in the hadronic environment. A repeated 
question on the nucleon decay raised in the past is how the extremely slow process of 
the lifetime of order 10^^ years or even larger is modified by strong interaction having 
the time scale of order, 10~^^ seconds. Thus one must deal with the short time limit 
of the decay law, and one must face the fact that the exponential decay law is not 
exact in quantum mechanics. To my knowledge even now there is no convincing 
calculation of the nucleon decay rate, assuming that the law of grand unified theory 
is known. 

In cosmology essentially every elementary process occurs in macroscopic cosmic 
medium. The process needs not involve a finite number of fundamental particles 
alone. For instance, one can consider a macroscopic quantum process such as the 
first order phase transition taking place in thermal environment. Here the order 
parameter of the symmetry, the magnitude of a Higgs field, is the system variable 
and it couples to essentially all light particles of the standard model that make up 
thermal environment. With these in mind, in cosmology there are plenty of places 
one has to estimate effects of medium. When one has a clear idea of how to separate 
the "small" system from the environment, the essential part of the problem becomes 
how elementary processes are modified in dissipative medium. We feel that the 
problem should be analyzed from the first principles of quantum mechanics, at least 
at the conceptual level. In condensed matter physics this problem is generically 
known as quantum Brownian motion or quantum dissipation. 

It has become increasingly clear to us that the bulk of the past works on quantum 
dissipation relies on the simplified form of quantum friction, namely dissipation which 
is local in time. This corresponds to the exponential decay law when one examines 
the fate of some initial excitation. This approximation is excellent in the most 



dominant phase of dissipation. But it fails both in early time and very late time 
behavior of the quantum system immersed in dissipative medium. Some formalism 
exists to deal with more general nonlocal dissipation, but they are not very useful 
to the problems we would hke to address. Indeed, some of our new results have not 
been recognized in the past. 

In this talk I shall first review what we have recently achieved with regard to 
the general framework of quantum dissipation. Then in the second part I shall 
discuss two applications of this formalism: particle decay process in thermal medium 
and parametric resonant effect in medium. These are related to some important 
problems in modern cosmology; baryogenesis, inflation, and the moduli problem in 
supergravity models. 

Most results described here are already in our papers, |I|], 0, [0. But I reorganize 
arguments with a new perspective, and add some new material. I attempted in this 
written version to explain basic results in a self-contained manner. 



II Basic Formalism of Quantum Dissipation 

There are two standard methods for the theory of quantum dissipation. The first 
one is the approach based on quantum Langevin equation |Q. The other is the path 
integral approach initiated by Feynman and Vernon p, [|], [^, [0. Both methods 
have merits and demerits, and one should use both in clarifying specific details of 
relevant problems. Before going to explain some fundamentals of these approaches, 
let me explain how one models the dissipation for a small quantum system, being 
caused by interaction with a larger environment. 

The basic idea taken by most approaches is that existence of a continuously 
many environment variables coupled to a finite number of subsystem variables is 
the essential part of dissipation. After all one does not know all details of the 
environment. Indeed, ignorance should be the privilege here. It is fortunate that 
dissipation, or at least something recognized as such, occurs once one makes no 
measurement of the environment. Thus one postulates that detailed modeling of 
the environment and its coupling to the small system should be unimportant to 
dissipative behavior of the system. After integrating out the environment variables, 
one should have only a few phenomenological parameters to describe dissipation. 



ultimately obtainable from experiments. Despite of this phenomenological nature 
one should base all auguments on rigorous quantum physics, and the modeling is 
inevitable. 

One is thus led to the theory of the simplest, yst the most fundamental model 
of quantum dissipation, the linear dissipation in open systems. In this approach 
one models the environment by a continuously many set of harmonic oscillators of 
some arbitrary spectrum and couples it to the subsystem via a bilinear term. Let 
the subsystem variable in question be denoted by q and the environment variable 
by Qa- For simplicity we assume that the subsystem has one degree of freedom, but 
it should be evident to extend it to any finite degrees of freedom. We often use a 
discrete label for the environment variable, but actually we assume that there are 
continuously many of them. The Lagrangian of our problem consists of three parts: 

L = L,[q]+LQ[Q]+U,,[q,Q]. (1) 

We take for the system-environment interaction the bilinear term: 

L, = \e- V{q) , (2) 

LQ = ^T.iQl-^lQl)^2L ^^(^'H-^'^'M)' (3) 

/•oo 

L-mt = -(lJ2 Ca^a -^ -Q duj c{uj)Q{uj) . (4) 

00c taken to be positive is the smallest of the environment frequency spectrum. 

As a technical aside, the renormalization of the frequency will be discussed later in 
appropriate places. Introduction of a counter term in relation to the renormalization 
as often discussed in the literature is actually the problem of how to relate the 
parameters of the theory to observable quantities. 

The strategy of subsequent investigation is to integrate out quantum dynamics 
of the environment, by averaging over its initial ensemble given by some mixed state. 

IIA Operator method 

The Heisenberg equation of motion for the subsystem and the environment vari- 
able is derived from the above Lagrangian, 

'<i+V'{q) = - / duo c{uj)Q{uo ,t) , (5) 

Q {u,t)+u;^Q{u,t) = -c{uj)q. (6) 



The environment part is first solved as a sum of tlie liomogeneous and tlie inliomo- 

geneous solutions: 

^/ X ^ / X / X r^ / X sm(ujt) 
Q{uj , t) = Qi{uj) cos{ujt) + Pi{uj) 



c(a;) /■* , . c{uj) rtf 



00 

I dr sin uj{t — t) q{T) — / dr sina;(r — t) g(r) 

2u; Jo ztu Jt 

--^ / rfrsmcu(t-r)g(r). (7) 

2uj Jo 

Here Qi{uj) ,Pi{uj) are initial environment operators at t = 0, and tf{> t) is some 
arbitrary time. The system evolution is then determined by the quantum Langevin 
equation P, 

'q+V\q) + 2 f dTaj{t-T)q{T) = FQ, (8) 

Jo 

FQ{t) = - r duo^)(hMe-'^' + h\{uo)e'^') , (9) 



lijJc 

i i"^ 



«/M = -o / dujr{uj)e-'^\ (10) 
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oo 



rfc.DM^!M_prf,z,(,)^!M. (11) 

luj J —oo 2iUJ 

bi{uj) and bl(uj) are the initial creation and the annihilation operators of the envi- 
ronment part of variables. D{uj) is the density of states per unit frequency. The 
fundamental quantity r{uj) that characterizes both the environment spectrum and 
its coupling to the system is called the response weight. The right hand side of 
eq.@ describes the random force by the environment, and it is characterized by the 
following correlation function, evaluated in the initial state of environment: 

/oo 
dujr{uj) cosu;(r-s) (2ni(t^) + l)cnv (12) 

Only the symmetric part of the correlation is written here. ni{uj) = bl{to)bi{uj) is 
the number operator of the initial environmental variable. For instance, in thermal 
state of temperature T = 1//?, 

(2ni(c^) + l)env = COth(^). (13) 

To go further, we shall take dynamics of a single harmonic oscillator for the small 

system. 



Even with this hmitation there are an important class of apphcations as we shall 
later discuss. It turns out that the entire potential of the system plus the envi- 
ronment can be diagonalized. Although it can be diagonalized via real orthogonal 
matrix, it is more convenient to attach some physically unobservable phase factors 
to this transformation. The diagonal oscillator coordinates, Q(li^), along with inverse 
relations, are then given by 

Q{uj) = Q{uj) - ^2ujr{Lo)F{uj - ^0+) [q- dJ -^ :^-^ Q{uj') , (15) 



P{u) = Pico) - ^2co r{u) Ficu - zO+) \p- du' ^| ^— — P{u') | , (16) 



, puj'r{u') 



oo 



OO 



q = - dujJ2iur{iu)F*{iu-iO+)Q{uj), (17) 



Ulc 



p=- dujJ2ujr{Lu)F*{LU-iO-^)P{Lu), (18) 



oo 



Q{uj) =Q{uj) + J2uj r{uj) / dJ ^ — — Q{uj') , (19) 



UJc 
OO 



^2uj'r{uj 


')F* 


ico'- 


zO+) 


UJ^- 


UJ'^ 


-zO+ 




^2uj'r{uj 


')F* 


tO+) 



P{uj) = P{u) + ^2uj r (u;) / du' ^ ^— P(^') . (20) 

The analytic function that appears here is 

This function has cuts along the real axis, lo > uJc and uj < —Uc- We assume 
the threshold Uc > such that there is a gap between the two cuts. The following 
discontinuity relation holds: 

F{u + iO+) - F{u - iO+) = i27ir{u)F{uj + iO^)F{uj - iO+) = i2nH{u) , (22) 

along this cut. An important property of F{z) is that in the cut 2;— plane this 
function is regular except on the real axis, and that it can be analytically continued 
into the other Riemann sheets by the discontinuity formula. In the second sheet 
F{z) can be shown to have simple poles. The location of these poles is given by 

z^ -ujl-2nG{z) + 27rir{z) = 0. (23) 

The analytic property thus derived is based on the presence of the gap at — cjc < 
uj < ujc- by analytic extention through this gap the boundary values at cj ± iO+ with 
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— oo < u; < —ujc and ojc < oj < oo become related. Without the gap as in the case of 
ujc = 0, there is no unique way to analytically extend G{z) (hence F{z)) with a given 
r(u;) at cj > 0. Indeed, the gapless case is difficult to deal with in complete generality. 
In the case of odd regular function for r{z) the analytic extention is possible. On 
the other hand, in cases with the spectrum weight behaving like r{uj) oc cu" (with 
a a positive non-integer) near uj = 0, the simple analytic extention in terms of the 
two relevant Riemann sheets is impossible. 

What happened to the subsystem spectrum? The equation that governs the real 
part of its energy eigenvalue squared A is given by 

/(A) = A, f{X)^u;^,+V Tdu^-^, (24) 

with V denoting the principal part of integration. The function /(A) is related to 
the analytic function by 

/(A)-A = 3fJF-^(yA + iO+), (25) 

with 3? denoting the real part. In the \ < lo\ region one readily derives from 
the monotonic behavior of /(A) the condition for absence of isolated spectrum at 

< CU < CUc , 

/(^c~)>^c, (26) 

with ~ indicating the limit from below. Consistency further requires cuq > lOc- When 
this condition (^) is obeyed, the system frequency loq moves, with increase of the 
coupling c[iS) , into the second Riemann sheet of the complex uj plane across the cut 
starting at cj = lo^- It finally settles at the zero of eq.(p3D. The imaginary part of 
this zero then gives the decay rate of any initial excitation of the subsystem, as will 
be made clear shortly. 

When the response weight r(uS) — » a constant as cu — i> oo, as often happens and 
in field theory models later discussed in particular, one has to subtract a term from 
integrals containing the response weight. This corresponds to renormalization of the 
bare frequency loq with the frequency shift. 



oo 



The renormalized frequency is given by 



bu:^ = -2 r dJ-^. (27) 



oo 



r{ijj) 



o;;^ = cuq + 5a;^ = cjg - 2 / duj-^-^ . {2t 



The eigenvalue equation written using uj^ is then modified to 

X-ul-XV r du ^;^^^ =0. (29) 

The condition of stabihty of the entire system requires that the smallest eigenvalue 
A > 0, giving tu^ > , or equivalently 

ujI>2 r duj"^, (30) 

if the integral on the right hand side is convergent. In the discussion that follows 
the original subsystem frequency ooq appears only through the renormalized ujr or 
the real part of the zero (pS]) . In this sense cjq is not a fundamental parameter when 
physical consequences are discussed. 

The important quantity for physical interpretation of the results that follow is the 
overlap probability between the true eigen operator Q{uj) with the original system 
variable q, which is calculable from eq. ([T5|) : 



l((g|<5(^)))r = 2ujr{uj)\F{uj-iO+)\'^ = 2ujH{uj) 

= ^^^(^) (31) 

(u;2_^2_n(cu))2 + (7rr(cu))2- ^ ' 

We decomposed the analytic G{uj — iO^) into the real and the imaginary parts, 

27iG{uj-iO^) =U{uj)+i7ir{uj), (32) 

n(a;) =V r duj' ^^^ . (33) 

J ~oo 00 — iO 

In the weak coupling limit of c{uj) -^ the overlap probability (0) has the Breit- 
Wigner form with a complex pole at the zero of eq. ([23|) . As will be more evident 
shortly, it physically means that any system excitation is damped with a decay rate, 

27 = —^ , 34 

with bj the real part of the pole of Fi^z). 

With complete diagonalization it should be possible to explicitly write the oper- 
ator solution. The solution is 

q(t) = - r doj ^2ujr{uj)F*{uj - iO+) Q{uj , t) , (35) 

Q{uj,t) = cos{ujt)Qi{uj) + ^-^^^^^Pi{uj). (36) 

CO 
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The initial values Qi , Pi are then rewritten in terms of the original variables, q , p 
and 



00 ^ , -, I ^ , -. ,-|- 



hM = \li^QM + -7^P.{^), h\. (37) 



After some straightforward calculation one finds that 

<i{t) = Ptg{t)+qig{t) 

- r dujJ^) (h*{uj,t)e~''''k{uj) + h{uj,t)e''^'bl{uj)) , {3t 

Juic 

P{t) = P^9{t)+(ii 9{t) 



where we introduced 

/■oo 

g{t) = 2 dujH{uj) sm{ujt), (40) 

J UJc 

hiuj,t) = f dTg{T)e-'-\ (41) 





t 

iwt 



k(uj,t) = / dujg(T)e-"^^ = g(t)e~"^' + iujh(uj,t). (42) 

Jo 



These solutions, eqs.(^), (PUJ), are fundamental to any calculation of correlators 
such as {q{T)q{s))env ■ Moreover, it completely solves the quantum Langevin equa- 
tion, (H), albeit for the limited case of system dynamics of the harmonic oscillator. 
The operator solution consists of two parts, one that depends on the system initial 
operator values, qi ,pi, and the other that depends on the initial environment values. 
Thus the memory effect of initial system dependence is governed by the function 

git). 

To see more on the memory effect, let us first ignore the random force caused by 
environment, by taking the ensemble average over the initial environment variables 
that are assumed to have no expectation value, {Qiiyj)) = (Pj(u;)) = 0. The 
expectation value of the system variable then follows [§ 

^ = -nHt){q)-cit){p), ^ = {P), (43) 

fi^(t) = q-^, C{t) = '-^^PLI . (44) 



^^^t) = 3^.-^ , Cit) = ^-l^^ 
g g -g"^ gg -g"^ 



The quantities that appear in this equation are the time dependent friction {C{t)) 
and the time dependent frequency squared (f2^(t)), incorporating environmental ef- 
fects. 
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Actually the ensemble average is not necessary. Using the explicit form of solu- 
tion, one gets after eliminating initial g^ ,pi dependence, 

^ + n\t)q + C{t) ^ = - n\t) f, - C{t) /,-/,, (45) 

/,= r du,f{u)[h*{uj,t)e-'^'h,{u) + {\i.c.)) , (46) 



U= rfcuV^(r(cu,t)e-*-*6,(cc;) + (h.c.)) . (47) 



This may be viewed as a local form of quantum Langevin equation: the operator 
acting on the system variable q is the local differential operator with local coefficient 
functions of time. Note that this form is not available in the general case of arbitrary 
potential. In the general case the integro-differential equation is available. 
Let us note that the function g{t) has the behavior of damped oscillation, 

g{t)r^-fAxi{Cjt)e-^\ (48) 

UJ 

if one takes the Breit-Wigner form for H{uj), or the pole approximation for F{z). 
When the coupling is weak, this pole approximation is excellent in many practical 
applications. In this case 

fi2 (^) = ^2 ^ ^2 ^ ^2 ^ C(t) = 27 . (49) 

Thus the effect of the environment is very simple in the pole approximation: aside 
from the random fluctuating environmental force, the system experiences the fre- 
quency renormalization and a constant friction: 

fj^W-f.nli,^,. (50) 

There are however physical phenomena that the pole approximation fails to de- 
scribe, as will be shown later. Let us simply state here that both early time and 
late time behavior of the system variables are not correctly described by the pole 
appproximation. For instance, the correct C{t) = 0[t^] as t ^ 0^ and at late times 
both of Q'^{t) and C{t) decrease with powers of time. These are not obeyed in the 
pole approximation. Indeed, the pole approximation violates the positivity of the 
density matrix at early times, as will be shown later. 

The correct final behavior at t ^ oo is the power law of time, as seen in the 
following way. Using the discontinuity formula, one may rewrite the u integration 
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for g{t), eq.(ig), containing H{uj) = {F{uj + iO+) - F{uj - iO+) ) /(27rz) along the 
real axis into the F{z) integration, the complex z running both slightly above and 
below the cuts. The factor sin{ujt) is replaced by Q=e~*'^* in this procedure. A half 
of this complex contour can be deformed into the second sheet, and one thereby 
encounters simples poles in the second sheet. The intregral for g{t) may then be 
expressed as the sum of the pole contribution (at z = zq with Qzq < for a single 
pole) in the second sheet and the contribution parallel to the imaginary axis passing 
through z = Uc, both in the first (I) and in the second (II) sheet 0: 



(7(r) = 53 fc-*«^»^) e^^°^ + 3 / dye-y^ iFi{u, + ty)-Fn{uj, + ty)) 

^ ^ TT JO 



, (51) 



with K~^ = zo — ttG (zq) + inr'^ZQ) . 

As seen from this formula, the pole contribution given by the first term describes 
the exponential decay as mentioned already. The rest of contribution gives the 
power law decay at very late times; oc t^"^^ 0. The power — « — 1 is related to 
the threshold behavior of the response weight, r{uj) oc (u; — cUc)" . This continuous 
integral also describes the correct early time behavior. 

Emergence of the power law term in the quantum Brownian motion has been 
noted in some specific models, but we find this behavior as a general property in 
the presence of the non-local dissipation. We would like to stress that the analytic 
structure of the cut uj plane is very important to derive the power law decay. For 
instance, if one takes the Ohmic form of the response weight given by |]I| 

■■"V) - , , „.f^r, , „, , . (52) 



(CU2 - fi2 + ^)2 + ^2^ 



2 



then one does not obtain the power law decay. The point is that this function, 
consisting of pole terms alone, does not have the branch point singularity needed for 
the power law behavior. 

Quadratic quantities averaged over the initial environment ensemble are com- 
puted as PI 

POO 

{q\t))= / duj{2n,{uj) + l),^^r{uj)\h{uj,t)\'' 

+ 9\t) {p1) + 9\t) {qi)+ 9it)m { p,q, + q,p, ) , (53) 

{p\t))= / duj{2ni{uj) + l),^,r{uj)\k{uj,t)\^ 

J UJc 

+ 9\t) {pI )+ 9" {t){qj ) + g{t) 9 it){p,q, + q,p, ) , (54) 

11 



1 r°^ 

- { q{t)p{t) + p{t)q{t) )= du{ 2ni{u) + 1 )e„v r{u) h{uj , t)k*{uj , t) 

+ g{t)9{t) {pl)+m g {t){qj) + {g\t)+g{t) 9 {t))\ {P^Q^ + m) -(55) 

Here we assumed both {hi{uo)) = {bl{iu)) = and {hi{uo)hj{uo')) = {h\{uo)}yj{uo')) = 0. 
Equivalent expressions are obtained using identities: 



dujcoth.{——)r{uo)\h{uj,t)\^ = 2 dr ds g{t — T)aR{T — s)g{t — s) , {5Q) 



,f3u - rt 

2 ' ^ ' ' ^ ' " Jo Jo 

du coth{(-^)r(u)\k{u,t)\'^ = 2 f dr r ds g(t - T)aR{T - s)g(t - s) , (57) 
2 Jo Jo 

anir) = duo { 2n,{u) + 1 ) -^-^ . (58) 

Let us turn to the late time behavior. At asymptotically late times the initial 
value dependence disappears, as mentioned above, leading to 



q{t) ^ - dujJ^(F*{uj-iO+)e'''^'bi{uj) + {h.c.)) , (59) 

pit) -^ i r diOio^fi^ [F*{LO-iO+)e-''^'bi{Lo)-{h.c.)) , (60) 



since h{uj , oo) = F{uj — iO^). The system evolution is thus governed by the initial 
distribution of the environment variables, bi{uj) , bi{uj), with the following probability 
functions; 



,2 ^,., ,((g|Q(cc^)))'^ 



riuj)\hiuj,oo)\' = Hiiu)= '^^^'-"^ '"^ , forg(t), (61) 

2a; 

rico)\kiLU,oo)\' = Lu'Hico) = '^\{{q\Qiiu)))\\ for p(t) . (62) 

Note that the overlap probability | ((g| Q{uj))) p has a universal character, being 
determined by general properties of the system and the environment alone, irrespec- 
tive of their particular initial states. The overlap probability between the system q 
variable and the diagonal Q{uj) respects the unitarity relation, 

_ roo 

duj\{{q\Q{uj)))\'^ = duj2ujH{u) = l. (63) 

Asymptotic values of averaged quadratic quantities are 

(g2(oo)) = / duj{2n,itu) + l),^^HiLu), (64) 

roo 

{p\oo)) = / du;{2n,{u;) + l),^,u^H{ij). (65) 
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Taking the thermal bath with {2ni{uj) + 1 )cnv = coth(^) , along with the pole 
approximation, this gives the well known results: 

{q\oo))^^, {p\oo))^T. (66) 

Let us discuss more closely the asymptotic form of the occupation number in 
thermal bath, 

, , 1 ,p^(oo) _ 2/ N \ 1 

Z UJ z 

I f°^ , ,, ,Puj^ ,_ , CJ^^ ^^, ^ 1 

1 JUJc 

The temperature dependent part of the occupation number is then 



/ rf^coth(^)(cu+— )i/H--. (67) 



When the pole term dominates, or equivalently one approximates H{uj) by the Breit- 
Wigner function, n^ has the factor of Boltzmann suppression, e~ ^^'^ at low temper- 
atures. But the pole approximation is not good at low temperatures. Indeed, let us 
examine a typical example by taking the form of r{uj) = c (cj — cj^)" , in the range of 
Uc < uj < Q {Q ^ Uc) . We assume < a < 1 and u ^ Max ( cUc , T ) . The result is 



n 



— r(« + l)e-'^'^=T"+^ (69) 



iU^ 



r is the Euler's gamma function. 

This shows that instead of the exponential suppression at low temperatures, 
what is left in medium after the decay has a power-law behavior of temperature 
dependence (oc T"+^). An implication of this behavior to the unstable particle decay 
will be discussed in the next section. It means that the remnant fraction in thermal 
medium does not suffer from the Boltzmann suppression factor at temperatures even 
much lower than the mass of unstable particle. 

Related time evolution is as follows. The corresponding behavior of g(t) com- 
puted from the continuous part of H{uj) integration gives 

,(t)^--r(a + l) ;^^ ^ r (70) 

One can estimate the transient time t* from the exponential period to the power 
period by equating the two formulas of git) in their respective ranges, to obtain 

'* " ^ ^" ( 2cr(af 1)7-^0 ' ^''^ 
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with 7 = — Qzq the decay rate. For a very small c the factor inside the logarithm 
becomes large (oc c~^"~^), and by the time t^ the initial population has decreased 
like 

It may thus be claimed that the power law behavior is difficult to observe. But we 
shall later discuss that this may not be so in cosmology. 

We shall refer to our original paper on detailed discussion of correlation func- 
tions. But let me mention the asymptotic form. For both ti and t2 in the asymptotic 
late time region, 

r°° 0UJ 

{qiti)qih)) ^ dujH{Lu) cos LU ih - 12) coth{'^), (73) 

{p{h)p{t2)) ^ / dujuj'^H{uj) cos uj{ti-t2) coth{--). (74) 

The other correlator vanishes: ( q(ti)p{t2) + p{t2)q(ti) ) ^ . 

Finally let us mention how some of our fundamental formulas are related to 
quantities in complete thermal equilibrium. The limiting values of the analytic 
function G{z) as z -^ real uj are related to the well-known real-time thermal Green's 
function |T0[, when one takes the initial environment in thermal state. For a single 
oscillator, 

G^{u) = G{u + le) , G^{uj) = G{u - le) , (75) 

^» + l ^G\uj)=i / rfre-^tr (p^T [g(r) g(0) ] ) . (76) 



-I3ui 



The operator method is thus very powerful. On the other hand, it is difficult to 
extract quantum statistical nature of the subsystem state in this approach. For this 
purpose, the path integral approach is useful, to which we shall now turn. 

IIB Path integral method 

The basic idea is that one is interested in the behavior of the g— system alone 
and traces out the environment variable altogether in the path integral formula. In 
the influence functional approach by Feynman and Vernon one directly deals with the 
probability instead of the amplitude. This way one can compute the reduced density 
matrix that describes the state of the small system incorporating effects of interacting 
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environment. We define the influence functional by convoluting witli tlie initial state 
of the environment. To do so we assume for technical reasons that initially we may 
take an environment state uncorrelated with the system. The influence functional is 
thus obtained after integrating out the environment variables: 

T\q{r),q\r)] ^ J VQ{t) J VQ\r) J dQ, J dQ[ j dQ f J dQ)5{Qf - Q)) 

■K ( g(r) , Q{t) ) K* ( g'(r) , Q'{r) ) p,{Q, , Q[) , (77) 

K{q{r),Q{T) ) = exp{iSo[Q] + ^S,^,[q,Q] ) , (78) 

So [Q] + 5i„t [q,Q]= f dr ( Lq [Q] + L,^, [q,Q]) . (79) 

Jo 

The influence functional is a functional of the entire path of the system g(r) and its 
conjugate path g'(r). 

P^{Q^ ,Q'i}=J2 ^n C(Q. ')MQ^) , (0 < Wn < I) (80) 

n 

is the initial density matrix of the environment, which can be any mixture of pure 
quantum states with probability Wn- What deserves to be stressed is that one does 
not observe the flnal state of the environment, hence integration with respect to the 
flnal values of Qf = Q'f is performed here. 

Once the influence functional is known, one may compute the transition proba- 
bility and any physical quantities of the g— system by convoluting dynamics of the 
system under study. For instance, the transition probability is given, with introduc- 
tion of the density matrix p^^^ , by 



dqf J dq'frfiqf)p^''\qf,q'f)MQ'f), (81) 

p(^)= / Vqir) fvq'ir) j dq, / rfg^^^gO ^.(g.) ^[g(r) , g'(r) ] e^^'M-^'I'^'l , (82) 



where ipi,/s are wave functions of the initial and the flnal g— states, and Sq[q] is the 
action of the g— system. 

The form of the influence functional is dictated by general principles such as prob- 
ability conservation and causality. Feynman and Vernon found a closed quadratic 
form consistent with these, 

nQ{r),q'{r)] = 



exp 



- J^ dr J^ ds ( i{r)aR{T - s)i{s) + i e(r)a,(^ " s)X{s) ) ] , (83) 
with e(r) = g(r) - g'(r) , X(r) = g(r) + g'(r) . (84) 
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Thus two real functions ai{T) are all we need to characterize the system-environment 
interaction. These are defined here in the range of r > . The fact that «« depends 
on the difference of time variables, r — s, is due to the assumed stationarity of the 
environment. The Feyman- Vernon formula is valid for general Lq[Q] and Lg[q], not 
limited to the harmonic oscillator model if the interaction Ljnt [q , Q] is bilinear. 

The correlation kernels appear in the influence functional as a form of the nonlocal 
interaction and they are the dissipation aj and the noise an. The dissipation kernel 
a J thus computed agrees with the one deflned in eq.(]T(]|) |T|. Let us now specialize 
to the case of the oscillator bath of temperature T = 1/(3, which is described for a 
single oscillator of frequency u by 



/ \ 1/2 

yn coth^pcu/zjy 

UJ 

■exp 



iQ' + Q") cosh{f3uj) - 2QQ' 



2 smh{f3u) 
The noise kernel, along with the dissipation kernel, are then given by 



^5) 



anir) = ^^ T dcu cothAr{uj) e-'^^ (86) 

Z J —CO Z 

j poo 

aj{T) = -- dur{u)e-'--. (87) 

Combined together, it gives the real-time thermal Green's function: 

a(r) = anir) + laiir) =Y.cltT{ppT[Q{uk,r) Q{ujk , 0) ] ) , (88) 

k 

a{uj) = r drair) e"^^ = ^ V c^ (- L . ^I[L^ S{uj^ - ool)] .(89) 

As noted already, these are given in terms of the response weight r{uj), and are 
governed by the analytic function G{z). 

For the system dynamics we further assume a single harmonic oscillator of fre- 
quency ujq. In the path integral approach integration over the sum variable X{t) is 
trivial in this case, since both the local part and the nonlocal action above are linear 
in this variable: 

\ Jl [i{r)X{r) - ujlar)X{r)) 

- f dr r ds ( e(r)afl(r - s)^{s) + i ^{T)ai{T - s)X{s) ) . (90) 

JO Jo 
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Thus result of the path integration of the system variable X{t) gives the classical 
integro-differential equation for ^(t): 



d'^ 



+ ool e(r) + 2 J' ds as) ai{s - r) = . (91) 



f dr r dsar)aR{T-s)as), (92) 

Jo Jo 



dr^ 
The end result of the ^ path integral then contains an integral of the form, 

rt 
Jo Jo 

using the classical solution ^(r) with specified boundary conditions, ^(0) = C,i , ^t) = 

In the local approximation often used the dissipation kernel is taken to have the 
form of 

aj{T) = Suj^SiT)+7]S'iT), (93) 

with (5c(j^ representing the frequency shift and the r] term the local friction. This 
choice enables one to solve the ^ equation (|9ll) by elementary means. On the other 
hand, the noise kernel is usually given by the response weight of the form, 

r(^) = -^/(^), (94) 

IT \l 

with f{x) some cutoff function and Q a high frequency cutoff. The cutoff is needed 
to tame the high frequency integral of ^^(t). The simplest cutoff function f{x) = 
6(1 — x) gives an approximate form of aj{t) with the friction rj and 

Suj^^--nn. (95) 

71 

At high temperatures this approximation reduces to the well known classical form 
of the fluctuation, 

aR{r) = ^6{r). (96) 

n 

At low temperatures, howerver, the use of the cutoff function only in the noise kernel 
while retaining the local form of the dissipation makes validity of this approximation 
dubious. Our path integral approach [^ does not make this local approximation. 



and instead uses exact solutions of the classical equation ([9T|). This equation has 
been used in other approaches 0, too. 
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The rest of deduction uses the Laplace transform, and we shall be brief, leaving 
technical details to our original paper . Solution of the integro-differential equation 
(pip is, using (yf(r) defined by eq.(PO[), given as 

./N ^ ajt-r) ( g{t-r)g{t) \ 

with the dot denoting derivative. Both g{t) and g{t) can be shown to satisfy the 
integro-differential equation of the form {x = g or g), 

—^ + uj^x + 2 I dTai(t-T)x(T) =0. (98) 

at'' JO 

The reduced density matrix of the quantum system at any time is obtained 
from the action written in terms of the boundary values, Sci{^f ,Xf ; ^j ,Xi) , by 
CO nvo luting with the initial density matrix of the thermal environment. This action 
is computed as 

^^ci = -jef-j^--W^.^f+'-Xfif-'-X,i,, (99) 

U = 2 f dr r ds z{t) orIt - s) z{s) , (100) 

JO JO 

V = 2 f dr r dsy{T)aR{r-s)y{s), (101) 

JO JO 

W= f dr r ds{y{r)z{s)+y{s)z{r))aR{r-s), (102) 

JO JO 

y{r) = ^^, (103) 

ait) 

z{r) = ^(t-r)-^(t-r)^, (104) 

ar) - -6^-e/f.(t-.)-^^^). (105) 

9{t) V 9{tj J 

The master equation for the reduced density matrix, which holds for any initial 
state of the system, may be derived [0 from this effective action: 

^1 = -2|^ + ^eXp-.C(t)e| + 2.e^Z^,,(t)p + 4D.,(t)e||, 

(106) 



where some new quantites are defined as 

0„it) ^ Uil^U^±ll^Ml^W-^-,w), (107) 
^ \ 99 -g^ 2g g9 -g^ 2 J 

D^,(t) = U~gW+ll^M^w. (108) 

9 9 -g^ 
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The same master equation was derived by Hu, Paz, and Zhang 0, and their deriva- 
tion was later simphfied by ref |T^. These works however do not give exact solutions 



like ours. Moreover, our derivation is different, in that we deduce the master equation 
from the exact solution. 

We shall discuss the question of the positivity, using the master equation. This 
argument is based on that evolution of a pure quantum state necessarily requires 
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(|t.p^)_^^<0. (109) 

If this inequality is violated, the initial pure state with tr p^ = 1 would evolve to 
a state with tr p^ > 1. But with tr p = 1, this means that the state at t ~ has 
a negative diagonal element of the density matrix. Using the general equation for 
^tr p^ derived from the master equation, and the initial behavior, 

Cit) ~ -^(^(3)(o)2_^(5)(o)), (110) 

Dppit) ~ -^^-ij[l]<0, (111) 

D,,p{t) ~ f/~jX[l], (112) 

one concludes that 

(A 

We introduced the following notation: 



^^trpM ~-J[l](Ag)2t<0. (113) 



I[f{u) ]^ Tdu coth{f^)r{uj)f{u;) . (114) 

In this sense the correct approach satisfies the positivity constraint at this level. A 
more direct and general proof of the positivity will be given elsewhere |1^ . 

Comparison to the pole approximation, or the local friction approximation, at 
t ~ would be instructive. Here as the initial behavior one has 

C{t) = 27 , Dpp{t) = 0[t], D,p{t) = 0[t'] . (115) 

Hence at t ~ 

— trp2 ~ 27trp2. (116) 

This clearly violates the positivity constraint 

19 



For further discussion we take as the initial state a product of thermal states, 
a system of temperature Tq = 1/j3q and an environment of temperature T = l/f3. 
We may take Tq = T when we apply to the decay process of excited level initially 
in thermal equilibrium. On the other hand, in the limit of Tg — ^ it describes the 
ground state of the system harmonic oscillator. 

After a series of straightforward Gaussian integration we find the reduced density 
matrix as a function of X and ^, of the form, 

p'^^\Xf,^f) = 2^exp[-AXj-Bef + ^CXf^f], (117) 

A=—, B = -(h-^), C = — , (118) 

/i=J[|Mc.,t)n + ^coth(^)(^2 + c.oV), (119) 

h = m h{uj , t)k*{uj , i) ] + ^ coth(^) g (g +oolg) , (120) 

h = A I M^ ' ^) I' ] + ^ coth(^) (i W.-g") . (121) 

cjq is a reference frequency taken as that of the initial system state, and equated 
here to the initial bare frequency. If one so desires, either the renormalized cjr or 
the pole uj may be taken as another choice. But we imagine the situation a small 
system was added to a large environment at some time, its mutual interaction being 
absent prior to the initial time. In this circumstance it is appropriate to take cuq 
as the reference frequency. Since dependence on the initial state dies away quickly 
as time passes, the choice of the initial reference is not crucial for determining the 
behavior of states at late times. Both of h{uj ,t) and k{uj ,t) are already defined in 
the preceeding subsection. The density matrix p^^^ from which any physical quantity 
at time t can be computed has explicitly been given by the discontinuity, H{uj) or 
r{uj). 

The diagonal part of the reduced density matrix with q = q', namely with ^ = 0, 
is 

p^^\Xf = 2q,0) =2\ - exp[-AAq^]. (122) 



The width of the Gaussian peak of ~ l/y/A is a measure of how the system behaves. 
For instance, if the width increases with time, the system is more excited than 
originally prepared, while if it decreases, it is more deexcited. 
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The basic quantities that appear in the reduced density matrix are related to ex- 
pectation values of the coordinate and the momentum operators at the same moment 

by 

(.^) - ^-lu (123) 

{P') = 2^+1^ = ^3, (124) 

{^(qp + pq)) = ^ = h- (125) 
Thus one may write the density matrix as 

/^^ r 1 ^2 ( {p^) {qp + pq? \ ,2 I Aqp + pq) y , 1 ..^^x 

The reduced density matrix is thus characterized by expectation values of quadratic 
operators, just as in the case of pure Gaussian system without the environmental 
effect. 

It is sometimes useful to transform the density matrix in the configuration space 
to the Wigner function /H^(a:,p), 

f^{x,p) ^ / rfep(^H2x,0e-^^«, (127) 



= y— exp[-4^x^- ^^ 4g ]• (^28) 

The Wigner function is expected to give the probability distribution in the phase 
space {x ,p) when the semi-classical picture is valid. Expectation value of the number 
operator, namely the occupation number, in terms of the reference frequency, equated 
to the pole location ui here, is calculated most easily from the Wigner function: 

^"^ " ^~20d^^2^~2^~0^ ^A ^1QA~2 

It consists, except the trivial |, of two terms, the term lj/ (16 A) from the Gaussian 
width of the diagonal density matrix element and the rest from the kinetic term — ^. 
This formula of course agrees with that of the previous derivation in the operator 
method. We shall discuss the rate of particle production in more detail when we 
turn to the specific application, namely the case of periodic frequency uj'^{t). 
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An important measure to explore the behavior of the quantum state under the 
environment action is the effective entropy of the subsystem under study. Even 
if the entire system of the subsystem plus the environment is in a pure quantum 
state, the subsystem can have a nonvanishing entropy, since one cannot measure the 
environment and traces out its freedom. The reduced density matrix derived in the 
path integral approach is especially suited to calculation of the subsystem entropy 
|]I[. We define as usual the entropy of the system as (henceforth we omit the index 
R for the reduced density matrix) 

^ = -trplnp, (130) 

where the trace operation is performed on the system variable q. A useful device to 
compute the logarithmic matrix is to convert it to a power series according to the 
formula like 

The arbitrary power tr p" is calculable as 

tr p" = y ( n ^ii ) P(9i + <i2,qi-q2)--- piqn + qi,qn- qi) 

(4^)n/2 



Vb + Va)^ - ( Vb - Va)^ 



(132) 



A more efficient way to compute the entropy is first to observe, with the presence 
of the power series expansion, independence of the factor C in the entropy formula. 



This justifies neglect of the C term (by setting C = 0) in the density matrix ( 117 ) 
for computation of the entropy, and only for this purpose. The next step is to 
identify this density matrix with C = as an equivalent harmonic oscillator system 
of frequency uj under a thermal bath of temperature 1//5: 



, ^. exp(-^coth(^)e?-^tanh(^)Xn . (133) 

\ TT coth(/3c^/2) ^^4 ^2^^^ 4 ^ 2 ^ ^ J ^ ' 

The equivalence is possible only for B > A and is established by the parameter 
relation, 

u = aVAB, tanh^ = ^/9- (134) 
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With this identification the equivalent temperature and the entropy are given by [|l[ 



with 



X 



f ^ 4vag 

ln((v^+v^)/(v^-v^))' 

S = ^^-ln(l-e-'^^) 
e/J-^ - 1 ^ ^ 



2^/1/3-/1 = v/4(g2)(p2) - {qp + pqy . (137) 



The entropy 5* is a monotonic function of the single variable x = JB/A . Its limiting 
values are 

S -^ -^(x-1) Ini^-^, asx^l+, (138) 

— i> In — , as X — i> 00 . (139) 

It is interesting to note that the entropy initially increases irrespective of the 
relation between the two temperatures of the subsystem and the environment, Tq 
and T. This can be proved by first writing down the time evolution equation for the 
quantity x^ = 4 {q^){p'^) — {qp + pq)"^ : 

dx 

— = - 2C{t) x" - 16Z}pp(t) ( g2 ) - 8D,p(t) {qp + pq). (140) 

As t ^ 0^, the dominant term in the right hand side is the Dpp{t) term. Thus at 

t~0 

(It 

— ~2J[l]t(g2)>0. (141) 

This proves the entropy increase. 

When A <^B holds, considerable simplification follows: 



UJ 



4VAB, (142) 

t ~ 2i3, (143) 

S ~ - In ^ + 1 . (144) 

2 4^ ^ ' 

When B = A,OT equivalently {p'^){q^) = {qp^q) , the density matrix is factorizable 
into a product of '>p{q) 4'*{q') , since tr p" = 1. This is a pure quantum state, having 
the entropy S = . 
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When the initial subsystem temperature Tq = l//9o is nonvanishing, the entropy 
is 

initially. For Tq ^ cuq, Sj ~ In ^ . It ends finally with 

5/= In?, (146) 

to 

for B ^ A. As will be made clear shortly, this is equal to In ? in the high temperature 
limit. This variation, from ln(To/ciJo) to ln(T/ci}), is a result one would naively expect. 
Existence of the equivalent harmonic oscillator basis for S > ^ actually implies 
much more than a convenient means of the entropy computation. We have found 
the basis of the states that diagonalize the density matrix at any instant of time: 
it is the Fock basis | n ) of the equivalent harmonic oscillator of frequency u with 
diagonal elements: 

p„ = 2smh(^— )e-^'^("+2) = -^^ = — —= = . (147) 

Since the parameters, A and B, are time dependent, the diagonalization of the 

density matrix is instantaneous. One thus thinks of the subsystem temperature 

varying with time as 

1 x{t) 

however with frequency also changing as 

To explore a more detailed nature of the final state, let us work out the frequency 
distribution of populated harmonic oscillator levels. Suppose that one asks the prob- 

2 2 

ability of finding the system in the unit energy interval around a; = ^ + ^ g^ . We 
get the energy distribution function fioj) [0] from the Wigner function: 



UJ 



27r •'"^^'-' ^ 2 2 

2 [a f lo .AA AC^ + uj^. . . ^ ^ .^. 

— V^exp - — ( — + . ^° ) Io{ — z), (151) 
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Here 



vr JO 



■ w COS 9 



is the modified Bessel function, tliat beliaves as 

I .... , 



Io{w) 



as w — >• oo 



1 (as w — > 0) 



(153) 



(154) 



^/ 2ttw 

tuo is again tlie reference frequency not to be confused with the original subsystem 
frequency. 

This distribution differs from the thermal one and has some interesting features. 

At tu ^> oo, it decreases as 



D 



A 



y/ 27TUJ0UJ Z V B 



^z e 



-wD 



1 

^0 



4^ 4C2 + ul 
1 



4A , 4C2 + a;2 



2_^ 

AujoB ^ B 



1/2 



(155) 
(156) 



assuming 2 7^ 0. The high frequency tail of this distribution decreases more rapidly 
than that of the Boltzmann distribution, e~^'^. On the other hand, at a; = 



/(O) 



lA 



UJq \ B 

The average energy of the distribution is computed as 



(157) 



{u) 



dujuj f{uj) 



luqB f AA 4C^ + ujj 
AA V ^0 4:^0 B 



(15^ 



This average satisfies the relation. 



{u) = uJo{{n) + -) 



(159) 



when u is replaced by ujq in eq. (|129D . 

The exceptional case of z = with an exact exponential form is actually im- 
portant, since it is realized in thermal equilibrium. This is possible only when 
C = , and V AB = ujq/A , giving 






T 



T. 



(160) 
(161) 
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The normalized spectrum shape written in terms of a scaled energy, x = uj/{uj), 
is characterized by a single parameter S, 



fo{x-S) = -e-^/''lo{^^-^x) (162) 



1 -x/52 wVl-52 



"Z"^' Trf^expf-a;^-— cos^l, (163) 

Jo 0^ 



with which the energy distribution is given by 

f{uj) = -^fJ^;6] , (164) 



A ( AA 4C2 + cj2 \ ^ _ Wo ^^{q^){p^)-{qp + pqY 



X 



^ ^^^^ icuo^ Au,B I {v^)+ujl{q^) 2(n) + r ^^^^^ 

This function has the following limiting behavior: 

fo{x]6) -» e^^'astJ -^ 1, (166) 

h{x]6) -^ ^^e-^/^as^ -^ 0. (167) 

\'2'KX 

The first moment x happens to be the same in the two limits of 5 = 1 and 6 = 0. 
The main difference of these two limiting distributions is in the low x part. We call 
S the spectral shape parameter. 

The asymptotic late time behavior of the reduced density matrix is determined by 
h{iJ , 00) = /q°^ dr g{T)e~'^'^'^ , which is equal to the boundary value, F{u — iO~^) , thus 
giving the already known relations, r(a;)|/i(a' , oo)p = H{uj) , r {uj)\k{(jj ,oo)\'^ = 
uj'^H{uj) . Thus, at asymptotic late times 

A ^ Ufdcocoth{^)H{co)) , (168) 



B 



Equivalent ly. 



1 f^ Bill 

- / ducoi\i{^)uj^H{uj), (169) 

C ^ 0. (170) 

(g2) ^ r duocoi\i{^)H{uj) (171) 

Ju}c 2 

[■00 R, ; 

(p2) ^ / ducoi\i{^)uj^H{u). (172) 

JuJc 2 

These agree with those in the operator approach. In these computations the analytic 
structure and only that is important. 
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In the high temperature hmit 

(g^)~^, {p')r^T. (173) 

The renormahzed frequency cj/j is defined in the preceeding subsection, eq.(|28D and 
is given by 

In this contribution the continuous integral during the period of the power law de- 
cay is numerically subdominant, suppressed by the factor 1/T^ relative to the one 
from the exponential period. On the other hand, at low temperatures the contri- 
bution from the threshold region, uj ^ Uc, cannot be ignored, giving the dominant 
contribution to the power law period, as noted in the preceeding subsection. 

It may be useful to write crude formulas for asymptotic values in the Breit-Wigner 
approximation, using 2uj H{uj) ^ 6{uj — uj) : 

(g2(oo))^^coth(^), (p2(oo))^|coth(^), ^(oo)^.;, (175) 



T(oo)^T, x(oo) ^ coth(^) , ^(oo)^ _r % . (176) 



,l3uj 2uuJo 

These are not precise, but they are nevertheless instructive to understand precise 
formulas for these quantities. 6{oo) < 1 for t^o 7^ Co. But if one takes the reference 
frequency equal to u, one gets 6{oo) = 1 , in complete agreement with thermal 
equilibrium. 

In closing this section, we would like to stress that the two approaches we de- 
scribed, the operator and the path integral, are both useful and complement each 
other. In both approaches the maximal use of the extended analyticity has been 
instrumental to completely clarify the details of the linear open system. The model 
of the linear open system is thus analytically solvable. 



Ill Particle Decay in Thermal Medium 

Let us apply these considerations to the decay of unstable particle; v^ — >■ X + X • 
We assume that the decay product x is a part of thermal components that make up 
the environment. The parent particle (f of mass M may or may not be in thermal 
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equilibrium with the rest of mediuni: since we focus on the late time behavior, the 
initial state dependence may be ignored. 

Let me first explain how the field theory may be incorporated in the framework 
so far discussed. The point is that the infinitely many environment variable Qa can 
actually be a composite field of more fundamental fields. We consider as the system 
variable qt a Fourier component of some scalar field ^{x) in a fully relativistic field 
theory. As a simplest field theory model that may describe the decay process, take 
a Yukawa type of interaction with a new bose field x(x); 

Linl = J d-^X Cint , C-mt = - - ^X^ , (177) 

with fi a parameter of mass dimension. V is the normalization volume. By identifying 
the Fourier-mode with 

qk e^'-' = ^^=i= ( a, + al, ) e^''' , (178) 

we deduce the environment variable Qk{^) coupled to it, 

I duc{iu)Q,{iu) = ^ 1^ d^xx\x)e~'^-^. (179) 

The continuous label uj of the environment variable Qk{^) is thus identified to the 
internal configuration of two body states with a given total momentum k. This 
field theory model was introduced in ref ||ll. Henceforth we regard the field x as a 
fundamental variable of the environment, assumed to be in thermal equilibrium. 

The correlation function akir), or more conveniently its Fourier transform ak{uj), 
can be computed as in the real-time formalism at finite temperatures. To lowest order 
in /i^ 



tak(T-s)=i(^f /■ ^s^Jk-s,^ ^^r.,2,^ ^^.,2, 



d'xe'^-Hi {t[x\x,t)x\0,s)]pp) (180) 



2' Jv 
is calculable in terms of the x^propagator in the momentum space, 

, r,' ,_^. - -I^ H-' ' P - ml) ]. (181) 

We omit the suffix x for the x~Kiass m^ in what follows. 

As in the T = field theory it is easiest to first compute the discontinuity in- 
stead of the full diagramatic contribution. Physically the most transparent way 
to compute the imaginary part of the self-energy diagram, 55n(ct;) , is to use the 
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analytically continued expression ||T5[ from the imaginary-time formalism. The ana- 
lytically continued imaginary part, or more precisely the discontinuity, directly gives 
the response weight via 

r(uj) = 2'^U(uj) = 2ujr(uj), Q U(uj) = — (U(uj - ie) - U(uj + ie) ) , (182) 

2i 

where T{uj) is interpreted as a decay rate in thermal medium. 



The result for the x loop diagram at finite temperatures is well known ||15|| , 
|]I|], and physically consists of two parts; the process </) <-h> x + X i^i the region of 
to > ^/k"^ + Am? and the other process '^ + X"^ X (forbidden when all particles are 
on the mass shell, but allowed in thermal environment) in < cj < fc, where m is the 



daughter x mass. The response weight r{uj) does not vanish for \uj\ > yk'^ + Am? 
and |ci;| < \k\ from the kinematics of the decay and the inverse decay of particles 
off the mass shell, with the constraint of the momentum conservation. Thus a gap 



of the spectrum exists in A; < |ci;| < y/c^ + Air? . The finite non-vanishing mass 
of the daughter particle (m 7^ 0) is important for the existence of the gap and for 
associated physical consequences that follow. 

Let us explain some details of the calculation of the response weight in the sub- 
threshold region oi \uj\ < k . One loop contribution from the imaginary-time formal- 
ism is given by |jl5| 



I? roo 

3n(a;) = -^ J dE{n{E)-niE + iu)) , (183) 



uj k Am? 

2^2V^":^ 



^^ = -±oJl--^^^ (184) 



where 



^(^) = ^^ (185) 



is the Planck distribution function of T = 1//3. Since 

n{E) - n{E + cu) = n{E){ I + n{E + to)) - n{E + uj){l + n{E) ) , (186) 



the imaginary part ( |183[ ) for \iu\ < k is a sum of the two contributions, X + V' ^ X 
and its inverse process that is allowed to occur in thermal medium. Note that ip can 
be off the mass shell: u"^ — k"^ ^ the y? mass^. The factor 1 + n represents the effect 
of stimulated boson emission. 
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The result of computation is now summarized. For ui > y/W + Arn? the response 
weight is M 



Note that r{uj) -^ a constant (= 2Mx decay rate in the rest frame of if), hence 
the frequency renormahzation or the subtracted form of u integral is necessary. For 
< uj < k only the second term in the bracket of Eq. (|187D contributes. 



A useful, and adequate approximation we exploit for subsequent estimate is the 
weak coupling scheme with correct threshold and asymptotic behaviors incorporated: 

1 . . 



F z) 



— z^ + cj^ — iTTr(z) 

In this approximation we replaced the u dependent real part of the self-energy 11 (cu) 
by a constant, hence by the constant pole location u. 

A quantity of physical interest is the fraction of remaining particles given by the 
occupation number at asymptotic late times, 

1 , P? 9 > 1 -Bi- cut 1 / N 

nk = -{^ + uj^ql )-^^ + ^- (189) 



for each k mode {uk = y k'^ + M^). The temperature dependent part of this quantity 

is 

"^1=1 duj^^;^-{uj, + -)H{u;,k). (190) 

Jo e'^'^ — 1 ujk 

One must sum over momentum k to obtain the number density of remnants. 

We shall limit our discussion here to the decay that occurs when the parent if 

becomes non-relativistic, 

^fc~M + |^»T. (191) 

This condition is relevant in interesting cosmological problems of the neutron decay 
at the time of nucleosynthesis and GUT X boson decay at baryogenesis |jl6| . 



Computation of the temperature dependent part of the occupation number nf 
may proceed by deforming the contour of u integration, in the same way as in the 
discussion of (?(r) in the preceeding section. There are then two types of contribution: 
one is the pole term that gives the usual Boltzmann suppressed contibution of e~ ^'^'= . 
When mode-summed, it gives the number density. 



(27r)3 ' 27r 

30 



This is the famihar Boltzmann suppressed formula. 

The second one is contribution from the continuous complex path that gives the 
power law behavior of temperature dependence. A part of this contribution in the 



region uo > y/W+^m? is analytically calculable by using r{uj) = -^ + 0[m'^] , valid 
for a small daughter mass m. It is 



n 



^ ^' r dtp r d- ^ 




2n^l67rM^ Jo A "^^ e^^ _ 


-1 


vr fi^T^ 




1440 M3 





(193) 

This calculation however ignores complicated logarithmic factors in r{uj) of eq.(|18^). 
We numerically computed all terms including the logarithmic factor in r{uj) 
along with 0[m^] corrections. It turns out that the total contribution is ten times 
larger than the analytic result above: in the m — * limit, 

n ^ 2A X IQ-^ t-— . (194) 

M3 ^ ^ 

The main part of this large contribution comes from |ci;| < /c. With a dimensionless 
constant introduced by /x = gM, this gives, relative to the photon number density 

7^-2xlO--(-A_)^^. (195) 



.r2 



We wrote here the numerical value using Gp the weak interaction constant of mass 
dimensions — 2 ( Gf^n ~ 10^^ ), as if it were relevant to the neutron decay. 

One may estimate the equal time temperature Teq at which the power contribu- 
tion becomes equal to the Boltzmann suppressed number density, to give 

^^— , ^^7xl0-i\ (196) 

M 33' T3 ' ^ ' 

taking as an example /i = 10^^ M , the weal interaction strength. This number is in 
an interesting range to affect nucleosynthesis, but we should keep in mind that we 
did not work out the relevant three body decay, n ^ p + e + v^ ■ 

The physical interpretation of the pole term is a conventional one in terms of the 
remnant created by the inverse decay X^ X ~^ Vi with all relevant particles on the 
mass shell, hence suppressed kinematically by the Boltzmann factor e"^'^^'^. On the 
other hand, the contribution from the continuous contour can only be interpreted 
as remnant particles far off the mass shell that may exist in thermal equilibrium. 
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The usual kinetic approach such as the Boltzmann-hke equation is based on the 
rates computed from S-matrix elements on the mass shell and gives the Boltzmann 
suppressed abundance in equilibrium for M ^ T. Our fully quantum mechanical 
approach yields a different result. 

We note that the local friction approximation is equivalent to the pole model 
(with identification of Uq + buP' = (JftzoY , r] = — 2Qzo ) that ignores the con- 
tinuum integral around the threshold. Hence the pole model, or the local friction 
approximation, fails to correctly describe the off-shell remnant. 

We shall mention another application of immediate interest in cosmology; the 



heavy X boson decay for GUT baryogenesis. It has been argued |T^ that there 
exists a severe mass bound of order, 

mx > 0[ax mpi] ^ 10^^ GeV , (197) 

to block the inverse process of the X boson decay so that generation of the baryon 
asymmetry proceeds with sufficient abundance of parent X particles. The usual 
estimate of the mass bound mentioned above is however based on the on-shell Boltz- 
mann equation. More appropriate formula in this estimate is our remnant number 
density, 

nx^O[2xlO-^]gl — . (198) 

mx 

(In a more realistic estimate one should consider the X boson decay into quarks 
and leptons. But for an order of magnitude estimate difference in statistics is not 
important.) With the GUT coupling of g'x/^T^ = 1/40, the equal temperature is 
roughly 

Te,^f. (199) 

Thus, already at temperature of about half of the X mass the Boltzmann suppressed 
formula is replaced by the power formula. The kinematical condition for baryogenesis 
must be reconsidered in view of our off-shell formula. 



IV Parametric Resonant Particle Production by 
Coherent Field Oscillation 

We now switch to an entirely different dynamical system under dissipative envi- 
ronment. Coherent field oscillation often appears in modern cosmology. Two main 
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examples worthy of explicit mention are inflaton oscillation that gives rise to the 
hot big band from the null state after inflation, and the moduli field oscillation that 
often appears in supergravity theories. 

We consiser the case in which a quantum system field if is coupled to a coherent 
field oscillation ^(t) oc cos{m,^t) . The coherent field may be viewed as an aggregate of 
zero-momentum particles with some kind of precise coherence, m^ is thus the mass of 
these bosons. Here this field oscillation is regarded as given, hence we do not discuss 
field damping due to particle production of the </?— field. Our main concern here is 
in the effect of thermal bath on the </)— particle production caused by ^—oscillation. 

Under this circumstance the harmonic oscillator variable is identifed to qk = 
, -^ ( ttk + alfc ) , where the creation and the annihilation operators are those 
of Fourier k— mode (oc e*^'^) of the system field fix). For the oscillator-system 
interaction we take the quartic coupling given by 

la'e^'^ (200) 

primarily because this case has been analyzed in considerable detail without tak- 
ing into account the environment effect. The /c— mode variable qk then obeys the 
evolution equation, 

( ^ + ^fc(^) ]lkiT) = 0, ujUt) = P + ml + g^, cos^im^r) . (201) 

^0 is the amplitude of oscillation. The standard form of this type of equation is called 
the Mathieu equation and is usually written in dimensionless units: 



—- + h + 26 cos{2z) qkiz)=0, 
az'^ J 






(202) 


P + ml 
z - m.T , h- ^ '^ +26, 
mi 


6 = 


g'eo 

Amy 


(203) 



It is well known that this quantum system exhibits instability in infinitely many 
band regions [|l^] of the parameters of {\k\ ,gC,o) , or {h,6). It has been realized 



that this instability gives rise to particle production |T^ . The most important recent 
development is that for large amplitude oscillation particle production and associated 
field decay is greatly expedited, typically ending in a time scale of O[10 — 100] 
times 1/m^ [l^, |jl9|, pO[, [^, |2^. However, in all initial investigations so far no 
systematic estimate of the environmental effect has been attempted. This is precisely 
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what we wish to do in the present section. Throughout this section we assume that 
the time scale of cosmic expansion is much larger than the oscillation time scale 
so that the cosmological expansion may be ignored. This is usually valid at the 
explosive stage of particle production. 

This dynamical problem belongs to the more general class of dynamics of time 
dependent harmonic oscillator. We shall briefly mention how much of the general 
formalism for the simple harmonic oscillator is modified in this case. There are 
two time dependences in this system: the nonlocal effect of the correlation a{T — s) 
introduced by the environmental interaction, and the time dependent frequency uj{t). 
The most important cases of applications can be analyzed using the approximate 
localized friction for ai{T — s), since we are primarily interested in the late time 
behavior. In this case one can explicitly perform the path integral in terms of a 
classical solution for qk{T). 

The notion of the localized friction arises by focusing on the late time behavior 
of ai{T). By the late time here we mean the asymptotic time that is much larger 
than any time scale of excitation in the environment; t ^ 1/ (maximal excitation 
frequency) ~ 1/^2. The limiting form is then given by 

aj{T) = - \ dujr{uj) sin(cur) -^ r]6'{T) + 600^ 6{t) , (204) 

Jo 

POD 

V = -2 rfrr«7(r) = 7rr'(0), for r(cu = 00) = , (205) 

Jo 

o^r[t) -^ - / c?cJCj/(— ) coth(— -) cos(cjr). (206) 

TX Jo \l 2 

The essence of this approximation is the low frequency truncation to the response 
weight, 

r(cj) =cjr'(0) = -cu. (207) 

TX 

Although the dissipation given by otiij^ is local by the time scale larger than 
l/f2 the cutoff scale, the damping time scale of the noise kernel olb{j) may differ. 
It can be shown that at high temperatures of T > fi the locality of the correlation 
holds for r > l/fi, giving 

aij(r)~^5(r). (208) 

TT 

On the other hand, at low temperatures of T < 1] the local approximation may be 
inadequate and a high frequency cutoff is needed. As already discussed, we replace in 
this case the response weight by multiplying some cutoff function ]{x\ for instance 
the simple frequency cutoff at c<j = fi with ]{x) = ^(1 — x), as is often practiced. 
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With aii^T) oc 5'{t), the exponent factor in the path integral becomes, including 
the mass renormalization effect, 

-lJldrX{r){^^+ujl{r))ar)+'-rijldrX{r)i{r) 

'xi\^-'-r^ ifXf sign(t) - Ij dr £ ds ^{T)aR{T - s)^{s) . (209) 

A remarkable feature of this formula is that a part of the action excluding the an 
term is a local integral described in terms of a renormalized frequency c<j/j(r) and 
the friction term oc rj. This much is enough to considerably simplify calculation of 
the transition amplitude of the g— system. 

Path integral over the system variable q is standard: as in the more general case 
the equation for the classical path of the difference path C,{t), which is obtained by 
functional differentiation with respect to X{t), is 

This equation is best analyzed by introducing a new y(r), 

y(r) ^ e^-/' ar) , (211) 

(|^ + ^M^)-j)l/M = 0. (212) 

The technique of the Laplace transform, which was very useful for the simple 
harmonic oscillator in the previous sectios, does not work here due to the time 
dependence C(;|.(r). This is another main reason we have to resort to the local friction 
approximation. 

Limiting our analysis to modes within the instability band of the modified Math- 
ieu equation, we note that two linearly independent solutions are either growing or 
decaying according to 

y(r) = e^™«" P{t) , or e' ^™«" R{t) , (213) 

where A > is the dimensionless growth rate and both P{t) and R{t) are periodic 
with period of 27T/m^. There are two independent solutions for y{t), u(t) and 
v{t) satisfying the boundary condition, ■u(O) = 0, v(t) = 0. These are linear 
combinations of the growing and the decaying solutions. With the normalization 
appropriately given, leading asymptotic behaviors of these are 

u{t) ~ e^™«* P{t) , v{t) = - !^M!^ ^ e- ^'"^^ R{t) , (214) 

u(t) 
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with both P{t) and R{t) bounded functions. 

Note that the growth rate A here is related to the solution of the modified Mathieu 
equation, eq.( |212| ), and not of the original one. Thus we should keep in mind that 
the growth rate A does depend on the friction rj: \ = \{ri) . Suppose for instance 
that the renormalized y9— mass vanishes: m? + 600"^ = 0. The rj term in eq. (pi2|) 
then has an effect of lowering the band level, since the h parameter in the standard 
Mathieu equation is modified to 

/i = /i--% = ^ + 2^--%, (215) 

Ami mi Ami 

in this case. Although it is in general difficult to estimate the influence of the friction 
f] on the growth rate A(?7), the general trend is obvious: the lowering effect of the 
band level tends to increase the growth rate. In the extreme case of a very large rj 
the h parameter becomes negative, and the parameter is clearly in the instability 
region. In perturbation theory we later work with, the leading A term is however 
independent of r]. 

The exponential growth of the original mode function gfe(r) is possible for 
2X{r])m,^ > 1]. We observe two competitive factors for the growth: the rate of the 
parametric amplification 2X{ri) m^ against the friction 77. As will be shown in subse- 
quent computations of physical quantities, the friction, when it is small, does act to 
diminish the parametric particle production, but does not wipe out the parametric 
effect. Thus, instead of the blocking factor, it is more appropriate to view the friction 
as the inverse time scale for the system to be driven towards thermalization. Hence 
if the friction is small enough, the parametric amplification never loses against the 
friction. 

In subsequent discussion we only present analytic computation of time evolution 
of physical quantities and leave results of numerical calculation to our paper [|l|. In 
analytic estimate we concentrate on the asymptotic late time behavior for which 
r ^ max. ( 1/r] , l/(Am^) ) . With the asymptotic behavior of the classical solution, 
u{t) and v{t), the three quantities that appear in the reduced density matrix ( |117D 
of the V9— system are 

i3 ^ --— , (217) 
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^ 2ir u{t) 2^2 M(t) 



The quantities K ,A ,C are defined by 

-ryt ^0 



K 

A,B,C) 
a{uj) 

b{u) 
c{u) 



B + e 



— cotiil 

2 ^ 2 ' 



1 ,^^^(0),2' 



uV2 



^(0) 



,(3uj. 



duj coth( — ) r(uj) {a ,b ,c) (u;) 
2 

u{t) 



as ---— e 
w(0) 



jttJS— 2(t_5) 



ds 
■u(t) 









(218) 

(219) 
(220) 
(221) 

(222) 
(223) 



They have asymptotic behaviors: 
A = 0[l], C<0[Max. (e 



-Xm^t g-W2^ 



i^ = 0[Max. (e-^^^^Se""*)]. 

(224) 
Thus, we conclude that A,B ,C behave asymptotically as 



^ = 0[Min. (l,e- 



- 2Xm^t+rit 



B = 0[l], C = 0[1]. 



(225) 



We omitted dimensional scales of these quantities such as ujq and f], since these differ 
from case to case, depending on the parameter range of T , ujq , m^ , t] . For the 
reference frequency we take uJq = cj^(O) : the initial value of the frequency. 

The formulas so far are valid for any response weight r{uj). We now specialize 
to the environment system that may be described by the approximate form of the 
localized friction, r{uj) = - uj x (high a;— cutoff). In particular, we are very much 
interested in the high and low temperature limits of various quantities. The high 
temperature here effectively means that T ^ Max. ( ujq , m^ ,r] ,Q) . 

Here we shall write only some typical results for the average energy, the entropy, 
and the distribution function: for asymptotic late times we find at high temperatures 



s 



uJo{n) 



2Xm/: 



' r/t 



^m(O)^ 1^ coy 2 u{t)' 



Inf^l^le-*/^ 



\m^ ■u(O) 



-c^/(2M) 



2n{uo) 



(226) 

(227) 
(228) 
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Both the number (n) and the average energy (uj) increase exponentially with time, 
the rate 2A?7ig — rj being somewhat diminished from the value taken without envi- 
ronment. Note also that a rather non-trivial enhancement factor oc T appears in 
the prefactor. 

Despite of the suppression of the growth rate from 2Xm^ to 2A?Ti^ — rj, the envi- 
ronment interaction does not erase the parametric effect if 

2X{'r])m^>r]. (229) 

Let us look into the meaning of this condition that guarantees the exponential rate 
of particle production in thermal bath. Since r] is the relaxation rate of y?— field 
disturbance towards thermalization, the condition simply implies that the thermal 
relaxation never catches up the parametric amplification for a sufficiently small fric- 
tion. At the same time the condition implies that there exists a critical strength 
of the friction rj above which the parametric effect does not occur. The equation 
2A(?7) m^ = 7] thus gives the critical strength of the friction constant rj^. We wish to 
stress that although the suppression of the growth rate might have been anticipated 



from a naive consideration [|I^, it was not clear in the past how one can verify the 
parametric amplification itself, starting from quantum mechanical principles. In our 
view the significance of the friction rj is in its role as the relaxation rather than 
destruction of the coherence. Hence if the relaxation time scale is larger than the 
amplification time scale, the parametric effect wins over the decoherence due to the 
environment. 

We carried out detailed numerical analysis in |l[]. It supports analytical estimates. 
We only quote our main conclusions from that work: 

• For the first time it was verified that the parametric amplification can occur in 
medium, using the influence functional method. The average energy and the 
total number of population exponentially increases with time, if the growth 
rate of the parametric amplification in medium is larger than the friction r/, 
interpreted to be a measure of the relaxation rate towards thermalization of 
the system variable. The critical friction for the exponential growth has also 
been given, above which the parametric effect does not occur. 

• The resulting energy distribution deviates from the thermal Boltzmann dis- 
tribution, having a long tail of the high energy component characterized by 
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an exponentially increasing average energy and an exponentially decreasing 6 
parameter (spectrum shape parameter). 

Late time behavior of the average energy, the entropy, and the energy spectrum 
is insensitive to the environment temperature if time is displaced in comparison 
at different temperatures. This is true, even including the T = vacuum. 
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